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Abstract 

This paper presents joint power allocation and interference 
mitigation techniques for the downlink of spread spectrum 
systems which employ multiple relays and the amplify 
and forward cooperation strategy. We propose a joint 
constrained optimization framework that considers the 
allocation of power levels across the relays subject to an 
individual power constraint and the design of linear re- 
ceivers for interference suppression. We derive constrained 
minimum mean-squared error (MMSE) expressions for 
the parameter vectors that determine the optimal power 
levels across the relays and the linear receivers. In order 
to solve the proposed optimization problem efficiently, we 
develop joint adaptive power allocation and interference 
suppression algorithms that can be implemented in a 
distributed fashion. The proposed stochastic gradient (SG) 
and recursive least squares (RLS) algorithms mitigate the 
interference by adjusting the power levels across the relays 
and estimating the parameters of the linear receiver. SG 
and RLS channel estimation algorithms are also derived 
to determine the coefficients of the channels across the 
base station, the relays and the destination terminal. The 
results of simulations show that the proposed techniques 
obtain significant gains in performance and capacity over 
non-cooperative systems and cooperative schemes with 
equal power allocation. 



1 Introduction 

The use of multiple collocated antennas enables the ex- 
ploitation of the spatial diversity in wireless channels, miti- 
gating the effects of fading and enhancing the performance 
of wireless communications systems. Unfortunately, due 
to size and cost it is often impractical to equip mobile ter- 
minals with multiple antennas. However, spatial diversity 
gains can be obtained when terminals with single antennas 
establish a distributed antenna array through cooperation 
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[T)-|[3]]. In a cooperative transmission system, terminals or 
users relay signals to each other in order to propagate re- 
dundant copies of the same signals to the destination user 
or terminal. To this end, the designer must employ a co- 
operation strategy such as amplify-and-forward (AF) [01 , 
decode-and-forward (DF) [3 4| and compress-and-forward 
(CF) 0. 

Recent contributions in the field have considered the 
problem of interference mitigation and resource allocation 
in the context of cooperative communications with relays 
ll6l- lfT8l . This problem is of paramount importance in 
wireless cooperative cellular, ad-hoc and sensor networks 
|[T0l[T5l that utilize spread spectrum systems. Prior work on 
cooperative multiuser spread spectrum DS-CDMA systems 
in interference channels has not received much attention 
and has focused on the assessment of the impact of multi- 
ple access interference (MAI) and intersymbol interference 
(ISI), the problem of partner selection [4. 8 1 and the bit er- 
ror rate (BER) and outage performance analyzes [9 1. Other 
related contributions on resource allocation investigated the 
capacity of ad hoc networks ATI , cooperative spatial multi- 
plexing |[T4l . power and rate control lTT6l[T7l , and schedul- 
ing |[T8"1 . There has been no attempt to jointly consider the 
problem of resource allocation and interference mitigation 
in cooperative multiuser spread spectrum systems so far. 

In this work, we study the downlink of spread spec- 
trum systems which employ multiple relays and the AF co- 
operation strategy. Specifically, we consider the problem 
of resource allocation and interference mitigation in mul- 
tiuser DS-CDMA with a general number of relays, which 
have been originally reported in [19|. In order to facili- 
tate the receiver design, we adopt linear multiuser receivers 
||2T1|241 which only require a training sequence and the tim- 
ing. More sophisticated receiver techniques are also pos- 
sible l2Tl [30l . We propose a joint constrained optimiza- 
tion framework that considers the allocation of power levels 
among the relays subject to an individual power constraint 
and the design of linear receivers. We derive constrained 
minimum mean-squared error (MMSE) expressions for the 
parameter vectors that determine the optimal power levels 
across the relays and the linear receivers. In order to solve 



the proposed optimization problem efficiently, we also de- 
velop joint adaptive power allocation and interference sup- 
pression algorithms. Specifically, we derive computation- 
ally efficient stochastic gradient (SG) and recursive least 
squares (RLS) algorithms that can be employed in a dis- 
tributed fashion. The proposed SG and RLS algorithms are 
employed to mitigate the effects of the MAI and the ISI, 
and to adjust the power levels, increasing the capacity of 
CDMA networks with cooperative diversity. These algo- 
rithms can be implemented in a distributed fashion, which 
means the mobile units compute the coefficients and em- 
ploy a low-rate feedback channel to update the coefficients 
for the power allocation. In addition, other SG and RLS 
algorithms are developed to estimate the parameters of the 
channels across the base station, the relays and the destina- 
tion terminal of the cooperative DS-CDMA system under 
consideration. The proposed algorithms are compared with 
non-cooperative and cooperative techniques without power 
allocation via computer simulations. 

This paper is organized as follows. Section 2 briefly de- 
scribes a cooperative DS-CDMA system and data model 
with multiple relays. Section 3 is devoted to the problem 
statement and the constrained linear MMSE design of the 
interference mitigation receiver and the power allocation. 
Section 4 is dedicated to the derivation of constrained adap- 
tive SG and RLS algorithms for the estimation of the pa- 
rameters of the receiver and the power allocation across the 
base station, the relays and the destination terminal. Sec- 
tion 5 is devoted to the development of adaptive channel es- 
timation algorithms for the cooperative system under con- 
sideration. Section 6 presents and discusses the simulation 
results and Section 7 draws the conclusions of this paper. 



2 System and Data Model 

Consider the downlink of a synchronous DS-CDMA sys- 
tem communicating over multipath channels with QPSK 
modulation, K users, N chips per symbol and L as the 
maximum number of propagation paths for each link. The 
synchronous DS-CDMA systems is considered for sim- 
plicity as it captures most of the effects of asynchronous 
systems with low delay spread Il24l [30l . The system is 
equipped with an AF protocol that allows communications 
in multiple hops using n r relays in a repetitive fashion. It 
should be remarked that other cooperation protocols such 
as DF can be employed without significant modifications, 
however, the AF has been adopted for simplicity and due 
to its lower complexity for implementation [3 |. We assume 
that the base station transmits data organized in packets 
with P symbols, there is enough training and control data 
to coordinate transmissions and cooperation, and the linear 
receivers at the terminals are perfectly synchronized. Since 
the focus of this work is on the resource allocation and in- 
terference mitigation, we assume perfect synchronization, 
however, this assumption can be relaxed in order to account 



for more realistic synchronization effects in the network 
The cooperative DS-CDMA system under consideration is 
illustrated in Fig. Q] 




Figure 1: Downlink of proposed cooperative multihop DS- 
CDMA system. 



The received signals are filtered by a matched filter, sam- 
pled at the chip rate and organized into Mxl vectors r b d [i] 
and r;, ri [i], which describe the signals received from the 
base station to the destination and from the base station to 
the relays, respectively, as follows 

K 

rbd[n] = ^ Qfc rf [n] D k h bd [n] b k [i] 

k=l 

+ VbdM + n bd[n], 



K 



Tb rj [m] = ^ a br j [m]D k h brj [m]b k [i] 



(1) 



k=l 



+ J 7fcr j H +n brj [m], 
rj = 1,2, ...,n r , i = 0,1, 



,P-1 



n = n r i + 1, m = n r i + r j + 1 

where M = N + L — 1, n(,d[i] and n br . [i] are zero mean 
complex Gaussian vectors with variance a 2 generated at 
the receiver of the destination and the relays, and the vec- 
tors r/ bd [i] and rj br [i] represent the intersymbol interfer- 
ence (ISI). 

The M x L matrix D k contains versions of the signature 
sequences of each user shifted down by one position at each 
column as illustrated by 



d k {N) 






4(1) 

dk(N) 



(2) 



where dk = [4(1), 4(2), d k (N)] stands for the 

signature sequence of user k, the L x 1 channel vectors from 
base station to destination, base station to relay, and relay 
to destination are h b d[n], h br . [n], h rj d[n], respectively. By 
collecting the data vectors in (Q]l (including the links from 
relays to destination) into a (n r + 1)M x 1 received vector 
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at the destination we obtain 



rbd[n] 




r ri d[m] 




. r r nr d[m] 





J2k= 



J2k=i a bdl n ] D khbd[n]b, 
d [m\D k h rid [m]b r k x 



m 

r\d r 



d [m]D k h rnrd [m]b r k nr [i] 



Rewriting the above signals in a compact form yields 



A' 



r[i] =Y J C k / H[i\B k \i]a k \i]+r } \ 



t\ + Hi 



(3) 



Fixing a k [i], computing the gradient terms of the La- 
grangian with respect to Wk[i] and equating them to zero 
yields 

w k , op t[i] = R^lilPeu^]' k = l,2,...,K, (6) 

where R[i] = £f =1 C k H\i]B k \i]a k \i]a» \i]B* \t\U h [*]Cf - 
(T 2 / is the covariance matrix and PcwH = ^KW r M] = 
Cfe?i[i]afc[i] is the cross-correlation vector. The quantities 
R[i] and PcwM depend on a k [i\. By fixing U7fc[i], 
computing the gradient terms of the Lagrangian with 
respect to afc, pt[i] an d equating them to zero, we obtain 
the following expression for the power allocation vector 



k=l 



where the (n r + 1)M x (n r + 1)L block diagonal ma- 
trix C k contains versions of the spreading sequences of 
each user, the (n r + \)L x (n r + 1) matrix H[i] con- 
tains the channel gains of the links between the base sta- 
tion and the destination, and the relays and the destina- 
tion. The (n r + 1) x (n r + 1) diagonal matrix B k [i] = 
diag(6fe[z] b r k ld [i] . . .b 7 k " d [i]) contains the symbols trans- 
mitted from the base station to the destination (&&[£]) and 
the n r symbols transmitted from the relays to the destina- 



afc.optW = (Ra k \i] + A/) 1 p ak [i], k=l, 2, 



K, 
(7) 



tion (Wj. 1 [i] . . . V k n [i]) on the main diagonal, the i 
1 vector a k [i] = [a^ d [n] a k d [m] 



-l)x 

[m\ . . . a™ d [?n\] T of the am- 
plitudes of the links, the (n r + 1)M x 1 vector r)[i] with 
the ISI terms and (n r + 1)M x 1 vector n[i] with the noise 
components at the destination. 



3 Problem Statement and Proposed 
MMSE Design 

We are interested in jointly designing a linear receiver and 
determining the optimal power levels across the relays sub- 
ject to an individual power constraint. Let us consider an 
MMSE approach for the design of the receiver for user k 
represented by a (n r + 1)M x 1 parameter vector w k [i] 
and for the computation of the (n r + 1) x 1 optimal power 
allocation vector a k [i] ■ This problem can be cast as 



arg min E 

w k [i],a k [i] 



[Wk^opt [i\,ak,opt ml 

subject to af [i]a k [i] = P A ,k, 



\b k [i\ - w% [i]r[t 
k = l, 2, 



where R ak \i] = £f =1 B» m H W%™k[i\w%\i]C k n\i]B k \i] 
is the (n r + 1) x (n r + 1) covariance matrix 
and the (n r + 1) x 1 cross-correlation vector is 

Pah [i\ = E[b k [i\B»\i}n\i] H cZw k [i\]. 

The expressions in (O and depend on each other and 
require the estimation of the channel matrix H[i], which 
is identical for each user as we are dealing with a down- 
link channel. The expressions in (|6]l and (Q require matrix 
inversions with cubic complexity ( 0(((n r + 1)M) 3 ) and 
0((n r + l) 3 ), should be iterated as they depend on each 
other and on user k. It should be remarked that the pro- 
posed optimization problem is non-convex and may present 
multiple solutions due to the joint optimization of parame- 
ters. However, the experience with the proposed algorithms 
suggest that the solutions may be identical because we did 
not notice problems with local minima or loss of perfor- 
mance under different initialization. A study of the op- 
timization problem is beyond the scope of this work but 
seems to be an interesting topic for future investigation. In 
what follows, we will develop algorithms for computing 
a k ,opt [i] , Wk, op t [i] and estimating the channel matrix 7i [i] . 



Proposed Constrained Estimation 
Algorithms for Receiver Design 
and Power Allocation 



The expressions for the parameter vectors Wk[i] and ak[i] 
can be obtained by transforming the above constrained op- 
timization problem into an unconstrained one with the help 
of the method of Lagrange multipliers [32] which leads to 

K 

A,. =E[\b k [i\-wf[i](J2Cin[i\B l [i\a l [{\+r l [i]+n 
i=i 

+ \(a?[i]a k [i]-P Al k), k = l,2, ...,K. 



In this section we present adaptive constrained SG and RLS 
~x estimation algorithms to estimate the parameters of the lin- 
ear receiver and the power allocation. A key feature of the 
proposed algorithms is that they can be employed in a dis- 
tributed fashion. The only information that needs to be sent 
via a feedback channel is the power allocation vector. 



4.1 



Adaptive Constrained Estimation and 
Power Allocation with SG Algorithms 



(5) 



In this subsection, we will develop SG algorithms for com- 
puting Wk[i] and a[i] recursively. Let us consider the pro- 
posed constrained optimization in resort to the method 



3 



of Lagrange multipliers [32] and express the following La- 
grangian 



(O into the constraint a k [i + l]a k [i + 1] = P A , k and then 
solve the following quadratic equation: 



C = E 



K 



\h[i}-wfm 5^CH[*]S,[i]a,[i] + f|[t]+n[ 



1=1 



a\ 2 + b\ + c = 0, 



where the coefficients of the equation are 



(13) 



(8) 



a = Pam, 
b = 2P A , k 



a\w"[i]C k n[i]B k [i]a k [i] 



where tbfc[i], H[i], and a[i] are parameter estimates of the 
receiver, the channel and the power allocation to be de- 
termined. Due to the nature of the problem, we need to 
jointly estimate these parameters. To this end, we will de- 
velop joint SG algorithms that can perform this task with 
low complexity. 

We consider the Lagrangian in (|8) and calculate the in- 
stantaneous gradient terms of it with respect to and 
a[i], respectively, as follows: 



af[t]Bni]#"[i]C£tMt]) 



c = a(af [i\B% [i]H H [i]C% w k [i]e[{\ + t»£ [i]C k H[i]B k [i]a k [i]) 

+ a 2 (tfaf \i]C k n\i]B k [t]B? [0i H W^w k W|e[i]| 2 ). 
The solutions of this quadratic equation are 



V£,,,.r 



K 



^Ci^fijBjHajH+fiH +«[*]) • M 



1=1 



K 



™" W E CiK[i] B i[i]a-t[i] + r,{i] + n[ 



i=i 



= -r[i e 



(9) 



vc K[l] = -b* \?{h H \i\c%w k \i][b k {i\ 
( K 

- w%[i] 22CiU[tiBi[%\ai[%\ + r)\i] + n\i] 



i=i 

H 



= -B^[0i [i]Cf w k {i]e[{\ + Aa fc [z] 



(10) 



where e[i] = b k [i] — w k [i]r[i] — b k [i] — 
*f[*1(EjliW[»]Bi[»]o,[i] + ff[»] + n[%\) is the 
error signal, which is a function of tbfe[i], H[i], and a[i]. 

Adaptive SG algorithms can be developed by using the 
expressions for the instantaneous gradients above and using 
them with SG descent rules ll32l . yielding 



w k [i + 1] = w k \i] - nVC&*[i] 
= w k [i] + fj,e*[i]r[i] 



(11) 



a k [i + 1] = a k [i] - aV£ & * r 



= o fc [i] + a(B%[i\H H {i\C%w k [i\e[i\ + Ao fc [t]) 

(12) 

where /i and a are the step sizes for the recursions for the 
receiver and the power allocation, respectively. Notice that 
in the recursion for computing a k [i], the designer needs to 
determine the value of the Lagrange multiplier A. There are 
two alternative approaches to that. The first is to substitute 



Ai = (-6 - b^b 2 - Aac) /2a, 
X 2 = (-b + b\jb 2 - Aac) /2a. 

j These solutions have to be computed for every time instant 
i and checked before substituting them into ©. 

The second approach to computing the power allocation 
and ensuring the constraint is as follows. The constraint 
is relaxed at first by making A = and performing the 
following recursion: 

a k [i + 1] = a k \i] + a(Bf \i]M B [z]Cf w k \i]e\i]). (14) 

This is followed by a procedure to enforce the individual 
power constraint for each user k, i.e., a k [i + l]a k [i + 1] = 
P A>k , which is described by 

f Aa fc [i] ^ 

a k [i + l]<r- P A , k a k [i + 1] (\/af [i + l]a fc [i + 1]) 

(15) 

The algorithms for recursive computation of w k [i] and 
a k [i] require estimates of the channel matrix W[i], which 
will also be developed in what follows. A comparison be- 
tween the two approaches for deriving SG algorithms will 
be illustrated via simulations. The complexity of the pro- 
posed algorithm is 0((n r + 1)M) for calculating w k [i] and 
0((n r + 1) 2 M) for obtaining a k [i]. 

4.2 Adaptive Constrained Estimation and 
Power Allocation with RLS Algorithms 

Specifically, we consider the problem of the previous sec- 
tion using an exponentially weighted least squares criterion 
and develop an RLS algorithm for the proposed task. Let 
us now consider the following proposed least squares (LS) 
optimization problem 



[tbfc[i] s ojk[t]] = aig min ^a' l \b k [l] - w% [i]r[l}\ 2 

w k [i],a k [i] |— ^ 

subject to a^[«]ajfe[«] = P A ,k, for fc = 1, 2, . . . , K, 

(16) 
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where a is a forgetting factor. The LS expressions for the 
parameter vectors w k [i] and a k [i] can be obtained in a sim- 
ilar way to the previous section and are given for each user 
by 

(17) 



W k [l\ 

a fe [j] = (^ afe [i]+A/)- 1 p a J«] 



(18) 



where R[i] = Y^i=i a r[l]r [I] is the estimate of the 
covariance matrix and Pc«H = S;=i a ' =l ^fcH r M i s 
the estimate of the cross-correlation vector, R ak [i] = 
Ej=i al=lB " Wi"[l\Ck w k [l]w%[l]C k H[l}B k [l] and 
PaM = T?i=i* l=i b k [l}BZ[l}H[l] H C%w k [l}. Thequan- 
tity A is the Lagrange multiplier and also plays the role of 
regularization term. Due to the difficulty of obtaining a 
closed form for this parameter, we will rely on a numerical 
solution for obtaining an appropriate value for it. The 
expressions in ( fTTI l and ( fT8l require matrix inversions with 
cubic complexity ( 0(((n r + l)M) 3 ) and 0((n r + l) 3 ), 
should be iterated as they depend on each other and still 
require channel estimates. 

Our goal now is to obtain a recursive solution to the ex- 
pressions in (fTTI l and cfT~8T > and reduce the required compu- 
tations. To this end, we will resort to the theory of adaptive 
algorithms [32] and derive a constrained joint iterative re- 
cursive least squares (RLS) algorithm. This algorithm will 
compute w k [i] and a k [i] and will exchange information be- 
tween the recursion for improved performance. In order to 
develop the algorithm, we fix a k [i] and compute the in- 
verse of R[i] using the matrix inversion lemma [32 1 to ob- 
tain tb/s[i]. If we define 3>[i] = R[i] then we can obtain the 
recursions 



k[i\ = 



L 3?[z]r[i] 



1 + a~ 1 r H [i]<&[i]r[i] 



(19) 



and 



= a-^[i - 1] - oT^k[i\r H [i]<f>[i - 1] (20) 

Using the LS expression in (fTTI l and the recursion p CH [i] = 
o^Pc-hV ~ !] + fe fcW r M we 8 et 

w k [i] = R [i]p c «[i] = a$[i]p CH [i-l\ + $[i}r[i]bl[i] 

(21) 

Using the expression in < f20b for <&[?], substituting in d2~TI > 
and manipulating the terms yields 

w k \i]=w k [i-l] + k\i]£*[i\ (22) 

where the a priori estimation error is given by 

S[t\ = b k [t\-wg[i-l]r[%\. (23) 

The derivation for the recursion that estimates the power 
allocation follows a similar approach to the computation of 
Wk[i). However, there are some difficulties related to the 
enforcement of the constraint and how to incorporate it into 



an efficient RLS algorithm. At this point, a modification 
is required in order to complete the derivation of the pro- 
posed RLS algorithm. This is because the LS expression 
in ( fT~8b incorporates a Lagrange multiplier (A) to ensure 
the individual power constraint, which is difficult to embed 
within the matrix inversion lemma. Our approach is to ob- 
tain the LS expression for the problem in (fT&t without the 
constraint and then ensure the constraint is incorporated via 
a subsequent normalization procedure. In order to develop 
the recursions for a k [i], we need to compute the inverse of 
R ak [i] . To this end, let us first define & ak = R ak [i] and 
employ the matrix inversion lemma [32] as follows: 



a- x * ah [^B»{i]4l H [t\C%w h [i] 



k 



i + a-^fMC^W^H^MBfH^Hcf^K 

(24) 

and 

* afc [z] = a- l * ah [i-l]-a- l k ak [i]w% {i)e k n{i]B k \i}^ ak [i-l} 

(25) 

Now a recursive equation for computing p ak [i] can be de- 
vised by relying on time averages as given by 

i> a M = aiW* - + [i\n H [{\e^w k [i\ (26) 

The substitution of ( f26b into the unconstrained LS expres- 
sion yields 



a k [i] = R ak \i]p ak \i] 

= «*„ fc [z]p afc [i -l] + b k [z]* Qfc \i]C^w k [i\ 

(27) 

Substituting (f25l l into the above expression and after some 
algebraic manipulations with the terms we obtain 

a k [i\ = ® ak [i - l]p ah [i - 1] 

- k ak [i]C k H[i]B k [f]# , [# afc [i - 1} + k ak b* k [t] 

(28) 

By further manipulating the above expressions we get the 
recursive equation for [i] 

a k [i]=a k [{\ + k ak [i\C >1 [ i l ( 29 ) 
where the a priori estimation error for this recursion is 

U [i] = b k [i] - af [t]U H \i]C?w k W (30) 

In order to ensure the individual power constraint 
a^[i]afe[i] = Pa k , we apply the following rule 



a k [i] <- P A , k a k \i](Ja k [i]a k [i]\ (31) 

The algorithms for recursive computation of w k [i] and 
a k [i] require estimates of the channel vector H[i], which 
will also be developed in what follows. The complexity of 
the proposed algorithm is 0(((n r + 1)M) 2 ) for calculating 
w k [i] and 0((n r + l) 2 ) for obtaining a k [i]. 
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5 Adaptive Channel Estimation Al- 
gorithms 

In this section, we describe adaptive SG and RLS channel 
estimation algorithms for the cooperative DS-CDMA sys- 
tem operating with the AF cooperation protocol considered 
in Section 2. The proposed algorithms are developed for 
use with the parameter estimation algorithms derived in the 
previous section for receiver design and power allocation. 

5.1 Adaptive SG Channel Estimation 

In this part we present an adaptive SG channel estimation 
algorithm for determining the parameters of the channels 
across the links comprising the base station, the relays and 
the destination terminal. In order to derive such channel es- 
timator, we first cast it as the following optimization prob- 
lem 



across the links comprising the base station, the relays and 
the destination terminal. In order to derive such channel es- 
timator, we first cast it as the following optimization prob- 
lem 



H]i] = argmin 



1=1 



a l - l \\r[l]-C k U\i]b k [l]a k [ 



(36) 



Due to the structure of the (n r + 1)L X (n r + 1) channel 
matrix H[i] and its relationship with the power allocation 
vector a k [l], we found that it is convenient for the deriva- 
tion and algorithm development to combine them into a 
(n r + 1)L x 1 channel estimate vector 



h[t\ = H\i]a k [l] 
and recast the optimization problem as 



H\i] = argmin _E[||r[il 
HU] 



b k [i]C k ~H[i]a k [i] 



(32) 



Mi] 



arg mm 
hh] 



In order to derive an SG channel estimation algorithm, we 
start with the description of a cost function associated with 
the optimization problem in (|32| > described by: 

C = E[\\r[i]-b k [i]C k ii[i]a k [i]\\ 2 ] 



i=i 



b k [l]C k h[i] 



(37) 



(38) 



The solution to the above optimization problem is given by 
k]i] = * hPh [z] (39) 



= E[(r[i\ -b k ]i]c k n]i]a k ]i]) H (r[q-b k [qe k n[{\a k ]i])^ 



where 



= E[r H m^-bl[^[i]H H [ l ]C^r[ l 

H 



l^i=i < 



- l b* k Wk 



H. 



(CfCfe)- 1 and p h [i\ 
It should be remarked that 



the matrix inversion in 3>h can pre-computed and stored at 
..//.. „ H ...«a.. H ~ the receiver for systems with repetitive spreading codes. In 

bk[i]r [i]C k U[i]a k [i] - a k [i]H [i]C k C k n]i]a k ]i]l rdeT tQ develop a recU rsive algorithm for estimating the 

(33) channel, we express p h [i] via the following recursion 



Computing the gradient terms of (l33l with respect to the 
(n r + \)L x (n r + 1) channel estimate matrix H[i], we get 



VC. 



-Cfr[i]of[i]6|EW+CfC fc «[t]a fc [i]af[i] 

(34) 



Using the above result on the gradient of the cost function 
and resorting to a SG optimization recursion, we obtain 



U\i- 



i] = •«[»] 

= U[q+v(C%C k M\i]a k \i]a%[i] 
-C%r\i}a%[i\b* k ]i]), 



(35) 



where v is a step size. This SG algorithm for channel es- 
timation works very well and can accurately determine the 
coefficients of the channels across the links comprising the 
base station, the relays and the destination terminal. The 
complexity of the proposed SG channel estimation algo- 
rithm is 0(((n r + I) ML)). 

5.2 Adaptive RLS Channel Estimation 

In this part we present an adaptive RLS channel estimation 
algorithm for determining the parameters of the channels 



p h {i\=a Ph [i-l]+b* k \i]efr[i\ (40) 
Substituting it into (f39b we obtain 

h[i] = ah[i-l] + bi[{\& h C%r[{\ (41) 

Once h[i] is computed, we need to apply a transformation 
in order to obtain H[i]. This is carried out by manipulating 
algebraically the relation in (l37t with the post multiplica- 
tion of a k [I], which yields 

h\i}a%[l}=U[i\a k [l}aZ[l] (42) 

Now computing the Moore-Penrose pseudo-inverse of 
a k [l]a k [I] we obtain the following relation 

ii\i]=h\i]a^[l](a k [l]a^[l}y 

where (■)' denotes Moore-Penrose pseudo-inverse 
This procedure for channel estimation also works well and 
can accurately determine the coefficients of the channels 
across the links comprising the base station, the relays and 
the destination terminal. The complexity of the proposed 
RLS channel estimation algorithm is 0(((n r + 1)ML)). 



(43) 
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6 Simulations 

We evaluate the bit error rate (BER) performance of the 
proposed joint power allocation and interference suppres- 
sion (JPAIS) algorithms and compare them with interfer- 
ence suppression schemes without cooperation (NCIS) [21 1 
and with cooperation (CIS) using an equal power allocation 
across the relays [8 1. We consider a stationary DS-CDMA 
network with randomly generated spreading codes with a 
processing gain N = 16. The block fading channels are 
generated considering a random power delay profile with 
gains taken from a complex Gaussian variable with unit 
variance and mean zero, L = 3 paths, and are normalized 
so that over the packets we have E[H H [i]H[i]] = 1. We 
adopt the AF cooperative strategy with repetitions and all 
the relays and the destination terminal are equipped with 
linear MMSE receivers.lt should be remarked that the noise 
amplification of the AF protocol is considered 0. The re- 
ceivers have either full knowledge of the channel and the 
noise variance or are adaptive and estimate all the required 
coefficients and the channels using the proposed SG and 
RLS algorithms with optimized parameters. We employ 
packets with 1500 QPSK symbols and average the curves 
over 1000 runs. For the adaptive receivers, we provide 
training sequences with 200 symbols placed at the pream- 
ble of the packets. After the training sequence, the adaptive 
receivers are switched to decision-directed mode. 



N=16, K=6 users, P=1500 




SNR 

Figure 2: BER performance versus SNR for the optimal 
linear MMSE detectors. Parameters: A = 0.02. 



In the first experiment, we consider the proposed joint 
power allocation and interference suppression (JPAIS) 
method with the MMSE expressions of © and ©. We 
compare the proposed scheme with a non-cooperative ap- 
proach (NCIS) and a cooperative scheme with equal power 
allocation (CIS) for n r — 1,2 relays. The results shown 
in Figs. |2] and [3] illustrate the performance improvement 
achieved by the proposed JPAIS scheme, which signifi- 
cantly outperforms the CIS and the NCIS techniques. As 



N=16, SNR = 12dB, P=1500 



10" 



O NCIS(n=2) 
□ CIS(n=1) 








* JPAIS(n=1) 




□ CIS(n =2) 




» JPAIS(n=2) 







2 4 6 8 10 12 

Number of users 



Figure 3: BER performance versus number of users for the 
optimal linear MMSE detectors. Parameters: A = 0.02. 



the number of relays is increased so is the performance, 
reflecting the exploitation of the spatial diversity. In the 
scenario studied, the proposed JPAIS approach can accom- 
modate up to 3 more users as compared to the CIS scheme 
and double the capacity as compared with the NCIS for the 
same BER performance. 



N=1 6, SNR = 1 2 dB, K=6 users, P=1 500 
10°^— r— . . . 







B JPAIS-IPC(n=1) 

* JPAIS-IPC(n=2) 











500 1000 1500 

Number of received symbols 



Figure 4: BER performance versus number of symbols. 
The curves for the adaptive SG algorithms using the so- 
lutions for the Lagrange multipliers are in solid lines, 
whereas those of the adaptive SG algorithms with a simple 
normalization are in dotted lines. Parameters: A = 0.02 
(for MMSE schemes), fx = 0.025, a = 0.015, v = 0.01 
(for adaptive schemes). 

The second experiment depicted in Fig. [4] shows the 
BER performance of the proposed adaptive SG algorithms 
(JPAIS) against the existing NCIS and CIS schemes with 
n r = 1 and n r = 2 relays. The aim of this experiment is to 
compare the proposed SG algorithms that utilize the solu- 
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tions for the Lagrange multipliers in the recursions ( ) with 
the simpler method that introduces a normalization of the 
power allocation vector (, and ). The techniques compared 
employ SG algorithms for estimation of the coefficients of 
the channel. From the results we notice that the SG recur- 
sions that obtain the values of the Lagrange multipliers via 
the solution of the quadratic equation have a slightly better 
performance than the normalization-based approach. This 
is basically due to a greater precision in the computation of 
the power allocation coefficients. For this reason, we will 
adopt this version for the remaining experiments. 



N=16, SNR = 12 dB, K=6 users, P=1 500 
10°? , 




500 1000 1500 

Number of received symbols 

Figure 5: BER performance versus number of symbols. 
The curves for the adaptive SG algorithms are in solid lines, 
whereas those of the optimal MMSE schemes are in dot- 
ted lines. Parameters: A = 0.02 (for MMSE schemes), 
H = 0.025, a = 0.015, v = 0.01 (for adaptive schemes). 

The third experiment depicted in Fig. [5] shows the 
BER performance of the proposed adaptive SG algorithms 
(JPAIS) against the existing NCIS and CIS schemes with 
n r — 2 relays. All techniques employ SG algorithms for 
estimation of the coefficients of the channel, the receiver 
filters and the power allocation for each user (JPAIS only). 
The complexity of the proposed algorithms is linear with 
the filter length of the receivers times the number of re- 
lays n r , whereas the optimal MMSE schemes require cu- 
bic complexity. From the results, we can verify that the 
proposed adaptive estimation algorithms converge to ap- 
proximately the same level of the MMSE schemes, which 
have full channel and noise variance knowledge. 

The fourth experiment depicted in Fig. [6]shows the BER 
performance of the proposed adaptive algorithms (JPAIS) 
against the existing NCIS and CIS schemes with n r = 2 
relays. The techniques NCIS and CIS employ RLS algo- 
rithms for estimation of the coefficients of the channel and 
the receiver. The proposed JPAIS scheme and RLS algo- 
rithms estimates the parameters of the channel, the receiver 
and the power allocation. The complexity of the proposed 
adaptive algorithms is quadratic with the filter length of the 



N=16, SNR = 12 dB, K=6 users, P=1 500 

10° » , 




500 1000 1500 

Number of received symbols 

Figure 6: BER performance versus number of symbols. 
The curves for the adaptive schemes are in solid lines, 
whereas those of the optimal MMSE schemes are in dot- 
ted lines. Parameters: a = 0.998. 



receivers and the number of relays n r , whereas the pro- 
posed optimal MMSE schemes require cubic complexity. 
From the results, we can verify that the proposed adaptive 
RLS estimation algorithms converge to approximately the 
same level of the MMSE schemes, which have full channel 
and noise variance knowledge. This indicates that the pro- 
posed RLS algorithms work very well and can determine 
the coefficients of the channels and the receivers. 



7 Conclusions 

This paper presented joint power allocation and interfer- 
ence mitigation techniques for the downlink of spread spec- 
trum systems which employ multiple relays and the AF co- 
operation strategy. A joint constrained optimization frame- 
work that considers the allocation of power levels across 
the relays subject to an individual power constraint and the 
design of linear receivers for interference suppression was 
presented. We then derived MMSE expressions and SG 
and RLS algorithms for determining the power allocation 
and the parameters of the receiver. We also developed SG 
and RLS channel estimation algorithms were also devel- 
oped to compute the coefficients of the channels across the 
base station, the relays and the destination terminal. The 
simulations showed that the proposed algorithms can ob- 
tain significant gains in performance and capacity over non- 
cooperative systems and cooperative schemes with equal 
power allocation. Future work will consider a study of the 
proposed algorithms and their extension to MIMO, OFDM 
systems, time-varying channels and limited feedback is- 
sues. 
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